Stability of the density-wave state of a dipolar condensate in a pancake trap 
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We study a dipolar boson-fermion mixture in a pancake geometry at absolute zero temperature, 
generalizing our previous work on the stability of polar condensates and the formation of a density- 
wave state in cylindrical traps. After examining the dependence of the polar condensate stability on 
the strength of the fermion-induced interaction, we determine the transition point from a ground- 
state Gaussian to a hexagonal density-wave state. We use a variational principle to analyze the 
stability properties of those density-wave state. 
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I. INTRODUCTION 

The recent realization of a Bosc-Einstcin conden- 
sate of chromium atoms 0, Q opens up the study of 
quantum-degenerate gases that interact via the long 
range, anisotropic magnetic dipole interaction. This is 
an anisotropic and long-range interaction that leads to 
the appearance of a wealth of new properties past those 
characteristic of systems with isotropic interactions Jj|. 
In the case of 52 Cr s-wave interactions are normally much 
stronger than the dipolar part, but it has been shown ex- 
perimentally that by applying a magnetic field the s-wave 
scattering length can be lowered to values comparable to 
the Bohr radius so that the dipole interaction domi- 
nates the system. In the regime where the dipole-dipole 
interaction is dominant the condensate is characterized 
by the existence of a metastable state whose properties 
depends on the trapping geometry [f| 0, 0,Jp, EH EH, 
as was experimentally demonstrated in Ref. [8(. A roton 
feature has also been predicted to exist in these systems 
for appropriate parameters. There would considerable 
interest indeed in accessing the associated roton instabil- 
ity [HI, ES E3] j as it is characterized by the spontaneous 
generation of a periodic density modulation in the con- 
densate [lj| . Unfortunately that state is unstable against 
collapse, but we found recently that this difficulty can 
be circumvented by the addition of a small fraction of 
non-interacting fermions, resulting in a significant stabi- 
lization of the bosonic system in cylindrical traps 
The simultaneous tr app ing of bosonic and fermionic iso- 
topes of chromium [171 ] also points to an interesting di- 
rections, with the possibility to observe novel ground and 
metastable phases in quantum degenerate polar boson- 
fermion mixtures. 

This paper extends our previous study of dipolar 
boson-fermion mixtures from cylindrical traps to pancake 
geometries, showing that dipolar bosons can be stabilized 



considerably by increasing the boson-fermion s-wave 
scattering length, and more importantly the density- 
wave states of dipolar bosons can be stable for typical, 
zero-field boson-fermion interaction strengths. Section II 
describes the induced potential created by fermions on 
bosons, and Sec. Ill discusses the general properties of 
the energy functional as a function of boson-fermion in- 
teraction using a variational Gaussian ansatz. Section IV 
presents an analysis of the transition point to the density- 
wave state for various combinations of trap aspect ratios 
and boson-boson contact interaction strengths. Section 
V introduces a density- wave ansatz in the form of a series 
of shifted Gaussians. The stability of these density-wave 
states is determined by varying the width of each Gaus- 
sian to seek minima in the interaction energy for various 
boson-boson s-wave interaction. We find a condition for 
the stability of these types of density modulated states. 
Finally, Section VI is a summary and outlook. 



II. FERMION-INDUCED INTERACTION IN 
POLAR CONDENSATES 

We consider a mixture of Nb dipolar bosons of mass 
mb and Nf single-component fermions of mass mj con- 
fined in a pancake-shaped trap characterized by a tight 
harmonic potential of frequency lo z along the z-axis and 
a softer harmonic potential of frequency luj_ in the trans- 
verse direction. A polarizing external electric or magnetic 
field, is taken to be along the z-axis . The dipole-dipole 
interaction between two bosonic particles separated by a 
distance r is then 



3z 2 \ 1 
Vdd{r) = fldd( 1 - — j -j , 
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where g,±d is the dipole-dipole interaction strength. In 
the mean-field approximation, the energy functional for 
the order parameter (f>(r) of the dipolar condensate can 
be expressed as 

gN b 



E 



(j)*(r)H (l)(r)d 3 r 



|0(r)| 4 ^r (2) 
r')\<j>(r')\ 2 d s rd 3 r'+E iad , 
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where 



H - V 2 

2mb 



[A 



2 <x 2 



y 



(3) 



is the sum of the kinetic energy and the trapping po- 
tential and A = lu±/lu z . The second term in the energy 
functional @ denotes the contact interaction between 
bosons, characterized by the strength g = Airfi^abb/rrib 
with being the s-wave scattering length, and the third 
term describes the nonlocal dipolc-dipolc interaction be- 
tween bosons. Finally, the last term Ei nd accounts for 
the fermion- induced interaction Vind (fc) between bosons, 
given in linear response theory [18[ as follows. 

We assume a contact boson-fcrmion interaction of 
strength gbf = 2irh 2 cibf /m r , where a,bf is the boson- 
fermion s-wave scattering length and m r = rribm / / (mb + 
to/) is the reduced mass. The boson- fermion inter- 
action energy has the form gbf J nf(k)n(—k)d?k, with 
n(fc) and nf(k) being the bosonic and fermion densi- 
ties in the momentum space. The linear response of the 
fermions to a bosonic density n(k) can be expressed as 
rif(k) = Vi n a(k)n(k), so that for the effect of the fermions 
on the bosonic energy functional is 



„ _ 1 9b f N b 
ind " 2 (2tt) 3 



Vi 



ld (k)n(k)n(~k)d 3 k. (4) 



The explicit form of the induced potential is 



Fermi sphere in momentum space, 
approximation, kf is given by [I7j ] 



k f £ z 



In this local density 



(7) 




where l z = \fh/{mb(jj z ) is the oscillator length in the z 
direction. 

Following Ref. [l6[ the induced potential Vi n d(fc) is 
given by 



9b f i 



V ind (k)={ 



where 



-1 



9bp 



n=l 

kfmf/{irhf 
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k > 2kf, 



three-dimensional 
k„ ~\~ k„, -\- 



is 



fermionic density of states, and fc 2 = 
the square of the fermionic momentum. The induced po- 
tential is attractive for very low momenta and goes to 
zero with increasing momenta. 

Throughout this paper, we consider the parameters 
corresponding to a 52 Cr- 53 Cr mixture, and fix the num- 
ber of fermions to be Nf = 10 3 unless otherwise stated. 
The zero-field bosonic s-wave interaction for 52 Cr is taken 
to be abb = 103ao, and a&/ = 70ao for the boson-fermion 
scattering length j2(| . 



Vi„d(fc) =9bfXf(k), 



(5) 



where Xf is the density response function It 
is related to the dynamical structure factor S(k,u>), 
which is the probability of exciting particle-hole pairs 
with momentum k out of the Fermi sea, by Xf(k) = 
— 2 du)'[S(k, uj')/uj'}. For a non-interacting single- 
component Fermi system we have 



(6) 



p<k f 
\p+k\>k f 



where p = \p\, kf is the Fermi momentum, the excitation 
energy is = pk cos 9/rrif + k 2 /(2rrif), and 9 is the 
relative angle between p and fc. This expression assumes 
that the fermions are locally free, so that there is a local 



III. STABILITY OF THE GAUSSIAN DIPOLAR 
CONDENSATE 

We now proceed to determine the stability of the dipo- 
lar condensate, using a variational ansatz in the parame- 
ter space of the fermion-induced interaction and dipolar 
strength. The variational wave function is taken as the 
Gaussian 



(r) = 



1 



: exp 



y 



2d 2 



z 

2d 2 



(9) 



with d, d z being the variational parameters and 
J d 3 r\<fi(r)\ 2 = 1. Substituting this Gaussian ansatz and 
its Fourier transform into Eqs. and yields the en- 
ergy Eg of the condensate as 



m b £ 2 z 



1 /l 



Eg= 2 i 2 +?? ' ' ~^ 
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F{r,- X ) + (J) {g < E<+9>E>+g md (E<-E>+E<)} 



(10) 



where rj = d z /d, and kf 



kfl z . The deriva- Ef , E> , E% , E^ , E$ and of the function F are straight- 



tion of algebraic forms of the interaction-energy terms 
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forward but lengthy, and their derivations are relegated 
to an Appendix. In Eq. (jTUJ) we have also introduced the 
effective three-dimensional dipolc-dipole interaction 



m b N b g d d 



53d 



The coefficients g < and g > find their origin in the 
momentum-independent contact interaction, given by 
the summation of the s-wave boson-boson scattering and 
the constant terms in the induced interaction They 
are given explicitly by 



9b s v 



47r i g dd 
9 



Finally 



5ind 



,1<1 



9bf v 



(11) 



487r5 dd fc| 



is due to the nonlocal induced interaction. We re- 
mark that the effective contact interaction consists of 
the boson-boson contact interaction as well as the 
momentum-independent part of the dipole-dipole inter- 
action and of the induced interaction. From Eqs. (fTTT 
and (jlip we have that for low momenta 



■ 9s= 3 



9 



(12) 



It is known that in boson-fermion mixtures without 
dipolar interaction , phase separation occurs when g s be- 
comes negative [22|, [5^, [24| . In this paper, in contrast, we 
restrict our considerations to the case g s > by chang- 
ing the boson-boson contact interaction g, so that phase 
separation does not take place. 

The energy functional (|I0[) was minimized with respect 
to r] and d z /£ z for various parameter values, with our re- 
sults summarized in Fig. [T] Without the fcrmion induced 
interaction, c/i n d = 0, the ground-state energy of the sys- 
tem is not bounded from below, and the strong trap- 
ping in z direction creates a local minimum in the energy 
landscape as a function of d z and rj [Fig.f (a)]. For fi- 
nite induced interactions <7; n d > 0, we find in contrast 
that the energy landscape is characterized by two min- 
ima, as shown in Fig. 1 (b), (c), and (d). In Fig. 1(b) the 
global minimum, which occurs for (j],d z /l z ) ~ (3,2.3), 
is a Gaussian state with narrow width in the trans- 
verse x-y plane. The additional local minimum close to 
(rj,d z /£ z ) ~ (0.2, 1.4) is a metastable state. 

With increasing boson-fermion interaction (?i n d, 
though, the ground-state energy corresponding to the 
global minimum at ?/ > 1 approaches that of the local 
minimum [Fig. 1 (c)], and these minima eventually reach 
equal energies at a critical value of g- m d . The new ground 
state past that point is characterized by the parameters 
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FIG. 1: Energy landscape E g obtained by the variational 
Gaussian ansatz Eq. as a function of g = d z /d and d~/l, 
at a constant dipolar interaction strength g^d — 30. In these 
figures the energy is cut off at E g = 1.5 from above and 
Eg = —0.5 from below for viewability. (a) Behavior of E g in 
the absence of the fermion-induced interaction gind = (b) 
Energy E g for g± n a = 0.02. The plot is characterized by pres- 
ence of two minima: (i) the true, narrow Gaussian ground 
state located at (g,d z /l z ) ~ (3,2.3); and (ii) a pancake- 
shaped metastable state located at (g,d z /£ z ) ~ (0.2, 1.4). (c) 
Energy landscape for <7i n d = 0.06. The energies of the ground- 
and metastable states approach each other, (d) The broad 
Gaussian metastable state at (g,d z /l z ) ~ (0.2, 1.4) eventually 
becomes the true ground state for gi n d = 0.07. In the figures 
the global (ground state) and local minima (metastable state) 
are indicated by a circle and cross, respectively. 



?y < 1 and d z ~ t z [Fig. 
Gaussian in x-y plane. 



1 (d)], that is, it is a wide 



Experimentally one may either vary the boson-boson 
s-wave scattering length while keeping the boson-fermion 
scattering length constant, or vary g^f with constant 
s-wave interaction g. In a boson-fermion mixture of 
chromium isotopes, typically g- m< \ ~ 10~ 3 is small, which 
corresponds to an energy landscape that resembles that 
of Fig. Q|b) and the broad Gaussian corresponds to a 
metastable state. Such a metastable state has been 
achieved experimentally in experiments by the Stuttgart 
group [|[ . A system that offers the potential to reach the 
regime of Fig. [ljd) is provided by a mixture of bosonic 
87 Rb and fermionic 40 K with the scattering length of 
rubidium atoms tuned close to zero. In that mixture, 
a zero- field scattering length a^j ~ 250ao [2l[ gives 
g ind - 0.2. 
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IV. TRANSITION TO A DENSITY WAVE 

The excitation spectrum of condensates dominated by 
a dipolar interaction is predicted to exhibit a roton min- 
imum EH- As a consequence, a bosonic density- 
modulated state in the x-y plane may arise as a local min- 
imum, and it may actually have a lower energy than the 



metastable Gaussian state and for high enough bosonic 
densities [lf| . This section discusses the transition to the 
appearance of such a state as a function of number of 
bosonic particles for various combinations of trap ratio. 

We proceed by introducing the new variational wave 
function that describes a density-wave structure with tri- 
angular symmetry, 



a 



cos 



nkox 
d z 



2 cos 



nkox 
~2d7 



cos n 



V3~k y \ 



2d z 



(13) 



where <f>(x,y,z) is defined in Eq. ([5]), n is an integer, and fco = kod z <C rj = d z /d. Substituting this trial wave 
function (fT3|) into Eq. (1), we find that the scaled excess energy of the density- modulated state relative to the 
Gaussian state 



z(k ,ai,a 2 ,a 3 ) 



2md\ 



is approximately given by 

e(fc , ai,a 2 , a 3 ) 



— (oj + 4a 2 + 9a 3 ) + — ^~ 



(Ed-W — Eg) 



(a\ + 2a ai + aia 2 + a 2 a 3 ) 2 V c g(k ) 



{a 2 + 2a ia2 ) 2 V cS (V3k ) + \{a\ + 2a\ + 4a a 2 + 2 ai a 3 ) 2 V cS (2fc ) 
2{a\d 2 2 + afaj + alal)V e g(V7k ) + (2a a 3 + al) 2 V cS (3fc ) 



(14) 



where = rj 2 £ z /d z The normalization condition reads a§ + 3(a 2 + a\ + a 2 )/2 = 1, and the effective transverse 
potential is 



V e $(k ) w < 



§ + ■9 < - \l TT k o crfcx (-?=)+ 5ind ( ) k 2 
2 



V2, 



fJ 



3 47rg dd 



£ z \ 2 (2k f d z y 



L-' 2 



fc < 2kfd z , 
fco > 2kfd z , 



(15) 



In evaluating Eq. (fT4|) we kept only the first four terms 
n = 0, . . . , 3 of Eq. (fT"3")) as the energy converges at the 
transition point. The magnitude of the error in that ap- 
proximate expression is estimated to be of the order of 
e = exp(— fcQ?7 2 /4). In the subsequent calculations we 
consider values of fco such that this error is less than or 
on the order of 1CP 6 . 

We numerically determine the variational parameters 
that minimize the excess energy e(fco, Oi, a 2l a 3 ) as a func- 
tion of g < for fixed and g 3 &- These results are summa- 
rized in Fig. [21 which shows the critical number of bosonic 
particles Nj, such that e has a minimum for fco =/= 0. By 
tuning the trap aspect ratio A = uj^/uj z to higher val- 
ues, the critical number of bosons Nj, required to achieve 
e < 0, that is, a transition from a Gaussian to a density- 
wave metastable state with lower energy, is lowered as a 
result of an increased n^. In addition, we note that the 



value of fco determined by the present variational calcula- 
tion matches closely the position of roton minima at the 
excitation spectrum of the condensate as in Ref . [HI . We 
also checked the energy with density wave having square 
symmetry and find that the triangular one always has 
lower energy. 



V. DENSITY WAVE STABILITY 

In the previous section we determined the transition 
to the density- wave instability, assuming the ansatz (|13[) . 
As already mentioned, this transition is related to the ex- 
istence of a roton minimum in the excitation spectrum. 
However, assuming such a density wave does not pro- 
vide an answer to the question of its stability against 
collapse. This is the issue that we address now, restrict- 
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0.26 



FIG. 2: Critical number of bosonic 52 Cr atoms needed for the 
transition from the Gaussian state to a hexagonal density- 
wave state for the trap aspect ratios A = 0.25 (solid curve) 
and A = 0.22 (dashed curve). Here t z = 0.15/xm. 



ing our considerations to the regime where the boson- 
fermion interaction strength is close to its zero-field value 
a-bf = 70ao. This corresponds in the case of 52 Cr- 53 Cr 
isotopes to gind ~ 10 -3 . This is the parameter regime 



characterized by a Gaussian metastable state, see Sec. II 
and Fig. []Tb). 

Once the transition to the density-wave state has oc- 
curred the condensate can also be described as a super- 
position of shifted Gaussian wave functions within the 
two-dimensional Z = x + iy plane. 



(j} dv ,(Z,z) 



1 



y/Mir 3 / 2 S 2 d : 



M 

J>xp 

.7 = 1 



{Z-Ztf 



2 1 



2f 2 



z 
2d 2 



where j is the index of the lattice site, Z 3 



(16) 
the lattice 



vectors generating periodic density modulations, £ the 
width of each density peak, M. the total number of den- 
sity peaks, and I the distance between neighboring den- 
sity peaks. It is given by / = 2?r/ko where ko is the 
position of roton minimum in the excitation spectrum of 
the condensate. 

Using the results of Sec. Ill to investigate the stabil- 
ity of each Gaussian in Eq. (|16p . we now show that the 
presence of fcrmions substantially stabilizes the density- 
wave state, noting that a stable density wave should be 
characterized by non-zero values of £,/d z and d z /£ z . 

From Eq. (fl"6|) , the total interaction energy is found to 
be [H 



m b £ 2 E int M 



h 2 



93d 

2 fd, 

3 U 



{9 



<E< 



9 > E? 



E<)}+f(l,£,d z ) 



(17) 



where we have assumed that £ -C I, i.e., that neigh- 
boring Gaussian peaks in Eq. (jTSJ) have little overlap. 
Here E< , E± , E^ , E^ , Ef and F have the same form as 
in the Appendix, but a is now a function of (^,d z ,6), 

a = J d 2 z cos 2 + t; 2 sin 2 8. The function f{l,£,d z ) is 
responsible for the particular geometry of the density- 
modulated state. Other terms, on the other hand, just 
arise from the energy of each individual Gaussian and 
we can thus apply the results of Sec. Ill to the present 
analysis. 

Assuming a triangular crystal and including only the 
interaction between nearest neighbors, we have 

f°° ~ ( k 2 f 2 \ ( kl\ ~ ~ 
/(J, = 6 J o V eS (k)expl j Jo I — J kdk, 

with k = kd z , the effective two-dimensional potential 
V e s(k) is defined in Eq. (| 1 5[) . and Jo is the zeroth-order 
Bessel function of the first kind. 

We now discuss the minimum of the energy functional 



£int(£, d z ) as a function of the effective strength of the 
contact interaction g Sl which is varied by changing g < , 
see Eq. (fT2")l . Without boson- fcrmion interaction, the in- 
teraction energy is a monotonically increasing function of 
£ with £mt(£ - * 0, d z ) — > — oo. As a result each gaussian 
of the density wave Eq. p^|) collapses. A typical exam- 
ple of the dependence of the interaction energy on £ is 
shown as the dotted curve in Fig. 3(a). With non-zero 
boson- fermion interaction energy, in contrast, the inter- 
action energy (|17|) exhibits a minimum at a finite £ as 
illustrated by the solid curve in Fig. 3(a) f or d z — 1.8£ z 
and g < = 0.17. The existence of that minimum implies 
the stability of each Gaussian, that is, the stability of the 
wave function ()16l) . 

We can determine the dependence of the gaussian 
widths £ on the s-wave boson-boson scattering length by 
changing g < . To do this we first minimize Eq. (fl7|) as a 
function of g < , treating d z and £ as variational parame- 
ters. Figure 3 (b) shows the resulting value £/£ z as a func- 
tion of g< for A = 0.25, N f = 1000 and l z = 0.15/jm. We 
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(a) f'z 
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(b) 9< 

FIG. 3: (a) Interaction energy £i n t of Eq. (117[l as a function 
of £/4 at A = 0.25, 4 = 0.15/im, d z = 1.54, and g < = 0.18. 
For gbf = the energy is unbounded from below, but for 
gbf = 70ao, a minimum appears at £/4 « 0.1. (b) Log-scale 
plot of the normalized width £/4 that minimizes Sint as a 
function of g < . 
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k A 

FIG. 4: Effective potential V c s as a function of k±d z for the 
parameters of Fig. |3Jb) for <; < = 0.2 (solid curve), and the 
critical value of g < » 0.14 (dashed curve). 



observe that £ decreases for smaller values of the contact 
interaction g < , and for g < < 0.14, the minimum-energy 
state corresponds to the collapsed state, £ = 0. Below 
that point the density wave is unstable. 

This behavior can be understood from the effective po- 
tential by hrst integrating out the z direction to obtain 
the effective potential (£ z /d z ) 3 V c e(k±d z ), where V e g is 
given by Eq. (TH]) and k\ = k 2 x + k 2 y . The width d z 
of the Gaussians is determined by the position of the 
minimum in Eq. (fT7|) . For g < < 0.14, V c g(k±d z ) is at- 



tractive for all high momenta, but at the critical value 
g < « 0.14, it approaches zero as kj_ — > oo, as shown by 
the dashed curve in Fig. [U Hence the energy is always 
minimized when the wave function has zero width in the 
x-y plane, i.e., £ — > 0. For g < > 0.14 the effective po- 
tential is repulsive both at low and high momenta, and 
attractive in-between as illustrated by the solid curve in 
Fig. 3J Subsequently the energy is minimized for a finite 
value of £. 

To find the range of parameter space characterized by 
the appearance of density-wave states, we need to con- 
sider, in addition to stability arguments, the transition 
point discussed in Section III. Tabic 1 summarizes the 
range of s-wave scattering lengths and critical numbers 
of 52 Cr atoms necessary to be inside stable density- wave 
regime. 



TABLE 1 


A 


4(/z) 


a>bb(ao) 


iV 6 (xl0 4 ) 


0.25 


0.15 


16-21 


> 0.65 


0.2 


16-19 


> 2.3 


0.25 


15-18 


> 5.0 



TABLE I: Tabulation of the scattering length of boson-boson 
contact interaction abb and critical number of bosons Nt for 
different values of aspect ratio A and oscillator length 4 inside 
the density- wave regime with Nf = 10 3 and gbf = 70ao for a 
mixture of chromium isotopes. 



VI. CONCLUSION 

In summary, we have analyzed the stability of a dipolar 
bosonic condensate mixed with non-interacting fermions 
in a pancake trap at T = 0. We found that the fermions 
help stabilize the condensate for a significant range of 
boson-boson and boson- fermion interaction strengths. 
We then investigated the transition of the system from 
a Gaussian-like to the density-wave ground state as a 
function of number of bosons, strength of the contact in- 
teraction, and trap aspect ratio. Our central result is the 
use of a variational ansatz to show that while in a purely 
bosonic system the density-wave state is always unsta- 
ble it can be stabilized by the admixture of even a small 
boson-fermion interaction. 

In particular, this study leads us to the conclusion that 
a pancake-shaped 87 Rb- 40 K mixture, which has a large 
boson-fermion s-wave scattering length, should be abso- 
lutely stable in the dipole dominated regime. By tuning 
the s-wave scattering length it is possible to reach a sit- 
uation characterized by the appearance of a roton insta- 
bility in the excitation spectrum, leading to the existence 
of a stable density-wave state. However, due to the small 
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dipolc moment of rubidium atoms the transition to this 
stable density-wave regime needs a substantial number 
of atoms, of the order of 10 6 . 

Future work will discuss the effect of the dipolar na- 
ture of the fermionic isotopes on the condensate - with 
chromium atoms in mind -, the existence and stability of 
density waves, as well as possible extensions to rotating 
systems. 
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APPENDIX A: DERIVATION OF THE ENERGY 
FUNCTIONAL OF EQ. (TlOj) 



In this appendix we derive the energy functional 
Eq. (fTTJ)) by using the Gaussian ansatz of Eq. ©. First 
we consider the dipolar interaction energy, 



J 



E, 



dd 



9ddN b 
3(2tt) 2 



dk 3 



fc 2 



1 I exp 



r 



{kidi + (ki + k 2 )d 2 } 



(Al) 



where k 



k 2 . Evaluating this integral gives 



the form of dipolar energy 



E 



dd 



with rj = d z /d and 



F(y) 



tan-^vV-l) 



(„2_ 1)3/2 -I)' 



(A2) 



Next we calculate the energy due to s-wave boson-boson 
interaction and fcrmion-induccd interaction. To achieve 



this goal we break the integral of the interaction energies 
expressed in spherical coordinates into two parts, E s = 



E< + E> and E, 



dd 



Ey d , according to 



oo f^kf 

dk = dk 
Jo 



dk. 



(A3) 



2k, 



For k < 2kf and in spherical coordinate, the interac- 
tion energy including the s-wave and induced interaction 
is given by 



J 



e: 



E:„j — 



ffdd^fc 

2tt 



2k, 



9bf v 



Ang 



,1,1 



E 

n=l 



_k_ 



2n 



1 



An 2 - 1 



exp 



-(d 2 z cos 2 8 + d 2 sin 2 



k 2 dk sin 6d6 



where g < is defined by Eq. (JTTJ) . This integral contain 
contribution from both the local and nonlocal part of 
the induced interaction. After integrating over k and 
rearranging the terms we obtain the expression 



e: 



with 

E<(a) 
Ef{o) 



E< d =g<E<+g ind (E<+E<), 



it/2 



48k] 



erf(v / 2fc/cr) 



sin 9d9, 



'2n 



</ 2 exp(~2k 2 a 2 ) 



sin 9d8, 



(A4) 



EH*) = 



~ 2" +1 / 2 1 
2^A, 4 n 2 -l 2fc 2 («-i) X 

r(n + 3/2) - T (n + 3/2, 2k)a 2 ) 

n 2n+Z 



tt/2 



sin 9d9, 
(A5) 



Here a = y d 2 cos 2 9 + d 2 sin 2 9 is a function of (9, d z ,d), 
crf(y), erfc(y), T, and 7 denote the error function, com- 
plementary error function, gamma function, and incom- 
plete gamma function, respectively. 
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Similarly, the contribution for k > 2k f to the interac- 
tion energy that contains both the s-wave boson-boson 



J 



interaction and the nonlocal induced interaction is 



2tt 



2k, JO 



r, 2 i, °° /r>; \2n 

g > _ 9b f v ( 2k f 



1 



^9dd ^— ' \ k J An 2 - I 

I 



cxp 



k 2 

-y(d z COS I 



d 2 sin 2 i 



k 2 dk smOdO, 



where g > defined in Eq. (fTTj) . In this equation the first 
term in the bracket stems from the s-wave boson-boson 
interaction and the last term from attractive nonlocal 
part of the induced interaction. Again after evaluating 
the integral over k we get 



E: 



E dd = 9 >E \ - ffind-Ef 



where 



(A6) and 



E>(a) 



erfc(V2fc /( 7) . 
— ^ — - sm 8d8, 



q °° 9k* 2 ( n+1) ri 2 



2ir 

n-2 

E 



An 2 - 1 



(— l) m 2 T " +1 (2fe / ) 2 



=o 2k f 



2n-2m-3 



(2n - 3)(2n - 5)....(2n - 2m - 3) 
I 



exp(-2fc 2 fC r 2 ) C r 2m 



sin #<i#, 



(A7) 



The total interaction energy is the sum of the contri- Eqs. (|A5|) . (|A7[) as the energy E g in Eq. (JTUJ) converges 
butions in Eq. |[A"2"]) . lfA4]) . (fM)) and is given in Eq. (JTOJ). with this inclusion. 
In this present paper wc took n = 1, . . . , 3 in the series in 
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